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ABSTRACT 

We study the correlation between the radio, mid-infrared and far-infrared properties for a sample 
of 28 blue compact dwarf (BCD) and low metallicity star-forming galaxies observed by Spitzer. We 
find that these sources extend the same far-infrared to radio correlation typical of star forming late 
type galaxies to lower luminosities. In BCDs, the 24/^m (or 22/^m) mid-infrared to radio correlation 
is similar to starburst galaxies, though there is somewhat larger dispersion in their q24 parameter 
compared to their qpiR. No strong correlations between the q parameter and galaxy metallicity 
or effective dust temperature have been detected, though there is a hint of decreasing q24 at low 
metallicities. The two lowest metallicity dwarfs in our sample, IZwl8 and SBS0335-052E, despite 
their similar chemical abundance, deviate from the average q24 ratio in opposite manners, displaying 
an apparent radio excess and dust excess respectively. 

Subject headings: galaxies: starburst — galaxies: dwarf — infrared: galaxies — radio continuum: 
galaxies 



1. INTRODUCTION 

More than 35 years ago, a correlation between the 
10/im mid-infrared (mid-IR) and 1.4 GHz lu mino s ities i n 
galactic nuclei was first noticed bv ivan der KruitI (|1971h . 
In 1983, the launch of the Infrared Astronomical Satellite 
(IRAS) brought a new era in collecting large samples of 
data in the infrared, and sparked numerous studies inves- 
tigating the above mentioned correlation in star forming 
galaxies. Both the radio and infrared emission arc closely 
related to star formation activities. The radio contin- 
uum is known to originate from two processes: thermal 
free-free emission from ionized gas in Hii regions and 
non-thermal synchrotron radiation from relativistic elec- 
trons that a re accelerated in the supernovae remnants 
(SNRs) (see iCondonI [1991 for a review). The ther- 
mal emission usually has a rather flat spectrum with 



f. 



oc f 



-0.1 



while the non-thermal component often has 



a much steeper spectral slope with f^ocz^""'^. In nor- 
mal galaxies the relative contribution of the two com- 
ponents varies with frequency and at 1.4 GHz the ra- 
dio continuum is dominated (at nearly ~90%) by the 
non-thermal component. The infrared emission is due to 
the thermal re-r a.diation of starlight from dust surround- 
ing Hii regions (jHarwit fc Pacini! I1975D . Based mostly 
on IRAS data, it has been shown that the far-infrared 
(FIR, 40-120 /im) to radio correlation holds well over a 
remarkably wide range of star forming galaxi es, spanning 



several orders of magnitude in luminosi ties (iHelou et al 
19851 Ide Jong et al.lll985l: iCondon et al.lll991l : lYun et al 



20011 ). The availability of deep observations of distant 
galaxies by the Infrared Space Observatory (ISO) and 
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recently by Spitzer (jWerner et al.|[200^ showed that the 
correlation between the IR and radio emission is not 
limited to the local univ erse but also extends in galax- 
ies a t higher redshifts (jGarrettl 120021 : lAppleton et all 
[200l . It was also revealed that in addition to the 
FIR, the mid-IR emission at ^24:^m correlates with 
the 20cm radio continuum, though with more scatter 
(Elb az et al.ll2002l:lGruppioni et~aLli2003l: lAppleton et all 
12004 IWu et al.ll2005f) . A number of very deep Spitzer 
mid-IR and FIR surveys can now probe a population of 
galaxies with low infrared luminosities for which ancil- 
lary data, inc luding deep radio iniaging, are becoming 
available (e.g. iJannuzi fc Devlfl999l: iSanders et ad] [20071 : 
[Rosenberg et al.ll20d6[ r This is particularly interesting 
since in the near future ~ 20/.tm is the longest wave- 
length which will likely be probed by the James Web 
Space Telescope. It is thus instructive to examine the 
mid-IR to radio correlation in more detail in these low 
luminosity nearby systems, for which little is known to 
date. 

As ubiquitous as the FIR/radio correlation appears 
to be, there are a few significant deviations from it. 
Radio excess exists in some instances, such as galax- 
ies hosting an active galactic nucleus (AGN). External 
magnetic field compression due to the interaction with 
nearby galaxies could also produce extra emission in 
the radio continuum (Miller & Owen 2001). Conversely, 
synchrotron deficiency has b een found in some nascent 
starburst galaxies studied by iRoussel et al.1 (|2003l l2006f) 
which was attributed to the lack of time for massive 
young stars to evolve into supernovae (SN) since these 
galaxies are just at the onset of a starburst episode. On 
rhi3@<ftftiie]!ik€a,ijiw®ii^cjaiti(ha)iiiia& 

damped in an optically-thin environment, because the 
ultraviolet (UV) and optical light may not be fully re- 
processed by the dust, as seen in low luminosity dwarf 
galaxies. 

Blue compact dwarf galaxies are a group of extra- 
galactic objects that are characterized by their blue op- 
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tical colors, small sizes (<lkpc) and low luminosities 
(Mb >-18). These galaxies do not display any AGN 
signature^ and have recent bursts of star formation in 
a relatively unevolved chemical environment. As such 
they have been proposed as nearby analogs of star for- 
mation in young galaxies in the early universe. In a 
metal poor environment, star forming regions are usu- 
ally optically thin and e mit le ss in the infrared. How- 
ever, iDevereux fc Ealej ()1989D suggested that in a low 
luminosity galaxy, the radio emission also decreases, and 
probably faster than the infrared. The deficiency in both 
the non-thermal radio and FIR emission may counterbal- 
ance each other and result in a similar FIR/ r adio ratio to 
the o ne observed in normal spiral galaxies (jKlein et al.l 
\imh . This was examined by a study of star forma- 
tion r ates (SFRs) in BCDs performed bv , Hopkins et al.l 
(|2002r ). in which the authors found an excellent agree- 
ment between the SFRs estimated from 1.4 GHz and 
60 /Ltm luminosities. As iBelll (|2003[ ) has pointed out 
though, the FIR/radio correlation is almost linear, not 
because the IR and radio emission reflect the SFRs cor- 
rectly, but because in low luminosity galaxies they are 
both unde restimated by si r nilar f actors. This is in agree- 
ment with iHelou fc Bicavl ()1993D . who found from their 
modelling work in disk galaxies that the transparency of 
the dis k was abou t the s ame to both re-emission pro- 
cesses. iHunt et"al] ()2005al ) in their analysis of the spec- 
tral energy distributions (SEDs) of low metallicity BGDs, 
noticed that these systems do not follow several of the 
usual correlations between the mid-IR, FIR and radio 
emission and display a scatter of a factor of '^10. 

The Spitzer Space telescope has enabled us to study 
the infrared properties of a large sample of BCDs, prob- 
ing the lower end of the luminosity an d metallicity 
range. In this paper, the fifth in a series (iHouck et al.l 
l2004bl : IW?rera!] [20061 . l2007a.b\ we examine their mid- 
IR and FIR to radio correlation extending the work of 
llfo^ ins et al.l ()2002f ). We describe the sample selection 
and the observational data in §2. A detailed study of 
mid-IR and FIR/radio correlation, as well as its depen- 
dence with other parameters, such as metallicity and dust 
temperature are presented in §3. We also discuss two 
extreme cases, IZwl8 and SBS0335-052 in §3. We sum- 
marize our conclusions in §4. 

2. OBSERVATIONS 

As part of the IRS^ (|Houck et al.ll2004al) Guaranteed 
Time Observation (GTO) program (PID: 85), we have 
compiled a sample of BCD candidates (~64) selected 
from the Sec ond Byurakan Survey (SBS), Bootes voi d 
galaxies (.Kirshner et all 119811: iPopescu fc Hopd |2000D . 
and other commonly studied BCDs. These sources are 
known to have low metallicities ranging from 0.03 Zq 
to 0.5 Zq'^. Their 22 fim fluxes have been published 

^ A possible exception among our sources is CG0752 in which 
the high ionization Unes [NcV] A14.3/24.3/xm arc detected (see 
[Hao et al. 2007). 

® The IRS was a collaborative venture between Cornell Univer- 
sity and Ball Aerospace Corporation funded by NASA through the 
Jet Propulsion Laboratory and the Ames Research Center. 

Here we ado pt the new oxygen so l ar abundance o f 
12-|-log(0/H)=8.69 fAllende Pri eto et al.l200l ). IWu etldl l|2007bl ) 
have derived neon and sulfur abundances for a sample of BCDs us- 
ing the infrared data, but these abundances are not available for 
all of the sources in this study. 



i n IWu et al.l (120061 ) . We also include 10 galaxies from 
Enge lbracht et al.l ()2005D (PID: 59), which are mostly 
BCDs and starburst galaxies, and span a larger metal- 
licity range (0.03 Zq'^1.5 Z©)^. For all galaxies with 
22 (24) /im detections, we searched the literature as well 
as the pubhc archives (NVSS and FIRST) for 1.4GHz 
radio continuum data. We restirct our sample to sources 
with both mid-IR and radio detections which results in 
a sample of 23 galaxies. Finally, we also include 5 galax- 
ies that have 22 ^m detections and 1.4 GHz upper limits 
published by Hopkins et al. (2002) for comparison be- 
tween ours and Hopkins' samples. Note that the sam- 
ple was not selected based on infrared properties, but 
merely on BCD-type objects and the availability of both 
infrared and radio data. As a result, our sample is not 
complete, but the large number of sources that only be- 
came detectable in the infrared with Spitzer, allows us 
to probe the properties of low luminosity dwarf galaxies, 
and provide statistically meaningful results. The obser- 
vational information for this sample and previously pub- 
lished data are presented in Table [U which includes the 
positions of the sources, their mid-IR, FIR and 1.4GHz 
flux densities, as well as oxygen abundances of the ion- 
ized gas. 

All sources in this study have mid-IR flux measure- 
ments either at 22/im with the IRS red peak-up camera 
and/or at 24 //m with MIPS (jRieke et al.ll2004D . The 
photometric fluxes of these two bands differ by less than 
10% for galaxies in the local universe. This was con- 
flrmed by using a suite of ^100 spectra of nearby galaxies 
from our IRS/GTO database and calculating their syn- 
thetic 22 and 24/im fluxes after convolving the spectra 
with the corresponding filter response curves. For con- 
sistency, in our analysis we use the MIPS 24 /xm mea- 
surements, and the 22 /xm values are only used when 
the 24 values are not available. For sources above 
the IRS 22/im and MIPS 24/im saturation limits we 
use the IRS low resolution spectrum to estimate a syn- 
thetic 24 /im flux. We also obtained far-infrared fluxes 
for our sample from the archiv al IRAS 60 and 100 /im 
data (Moshir & et al. 1990; SandersiranilQOl). 

Most of the 1.4 GHz radio contin uum data are from th e 
NRAO VLA Sky Survey (NVSS) (iCondon et al.l[T99l . 
while one is from the F ^aint Images of the Radio Sky 
at Twenty cm (FIRST) ([Becker et al.lll995[) . along with 
some individual observations (references in Table [J). A 
total of seven sources were too faint and were not in- 
cluded intheJA^EScataloCT For those we used the val- 
ues of iHopkins et al.l (|2002D who studied a similar BCD 
sample and remeasured the 1.4GHz fluxes using the im- 
ages from NVSS and FIRST, providing radio detections 
for 2 sources and better upper limits for another five 
sources that overlaps with the galaxies in our sample. 

Our final sample consists of 28 galaxies, all of which 
have Spitzer mid-IR 24/im and/or 22/im flux measure- 
ments. Among these galaxies, 23 sources have 1.4GHz 
radio continuum data and 5 have measured upper lim- 
its. IRAS 60 /im and 100 //m fluxes are available for 16 
sources and 3 more are detected only at 60 /im. We list 

* The galaxies in PID 59 partly overlap with the BCDs in PID 
85. We have also excluded some galaxies which the authors derived 
their 24 /im fluxes by doing a color correction to the IRAS 25 fira 
fluxes because of the aperture difference between these two band 
fllters. 
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TABLE 1 

Properties of the Sample 



ID 


Object Name 


RA 
J2000 


Dec 
J2000 




— 

Flux (m 




Jy) 




-I r* 1 1 / / T T \ 

12+log(0/H) 




References 




1.4 GHz 


24(22) /im 


60 fim 


100 fj.m 


1.4 GHz 


24 fim 


IRAS 


Z 


I 


Harol 1 




-33rl33m1 Qs 


26.8 


1900 


6880 


5040 


7.9 




(2) 


(3) 


(4) 


2 


NGC1140 


02h54m33.6s 


-10d01m40s 


23.6 


316.6 


3358 


4922 


8.5 


(1) 


(5) 


(6) 


(7) 


3 


SBS0335-052E 


03h37m44.0s 


-05d02m40s 


0.46 


66 






7.3 


(8) 


(2) 


(9) 


4 




04h30m47.0s 


+64d50m59s 


336.3 


2991 


54360 


55290 


8.2 


(1) 


(2) 


(3) 




5 




05h55Tn42 6s 


W O Li OHIO ^ o 


32.5 


1500 


6570 


5270 


8.1 




(2) 


(3) 


(7\ 


g 


UGC4274 


08hl3ml3.0s 


+45d59m39s 


10.5 


240 


3236 


6448 


8.5 


(1) 


(2) 


(6) 


7 


He2-10 


08h36ml5.2s 


-26d24m34s 


83.8 


4900 






8.9 


(1) 


(2) 


(11) 


8 


NGC2782 


09hl4m05.1s 


+40d06m49s 


124.5 


960 


9170 


13760 


8.8 


(1) 


(2) 


(3) 


02) 


9 


NGG2qn3 


09h32ml0.1s 


+21d30m03s 


444.5 


2200 


60540 


1 30430 


9.3 


(1) 


(2) 


(3) 


(11) 


10 


IZwl8 


09h34m02.0s 


+55dl4m28s 


1.83 


5.5 






7.2 




(2) 


(9) 


11 


RBS0q4n+544 


09h44ml6.7s 


+54dllm33s 


<2.3 


2.3 






7.5 


(14) 


(5) 




(15) 


12 


NGC3077 


10h03ml9.1s 


+68d44m02s 


29.0 


1500 


15900 


26530 


8.6 


(1) 


(2) 


(3) 


(16) 


13 


Mrkl53 


10h49m05.0s 


+52d20m08s 


4.0 


29 


280 


<480 


7.8 


(1) 


(2) 


(6) 


07) 

y^' J 


14 


VTT7,w403 


1 1 h27Tn5Q Qs 


-l-7Sd5Qm3qs 


1.2 


28 






7.7 


CIS"! 


(2) 






15 


Mrkl450 


Ilh38m35.6s 


+57d52m27s 


<2.0 


48 


279 


<575 


8.0 


(14) 


(2) 


(6) 


(7) 

y ' } 


16 


UM448 


Ilh42ml2.4s 


+00d20m03s 


32.6 


560 


4139 


4321 


8.0 


(1) 


(2) 


(6) 


(7) 

y ' } 


17 


UM461 


Ilh51m33.3s 


-02d22m22s 


<2.6 


30 






7.8 


(14) 


(2) 


(9) 


18 


UM462 


Ilh52m37.3s 


-02d28ml0s 


5.8 


110 


944 


896 


8.0 


(1) 


(2) 


(6) 


(9) 


19 


SBS1159+545 


12h02m02.4s 


+54dl5m50s 


<2.3 


6.4 






7.5 


(14) 


(5) 


(19) 


20 


NGC4194 


12hl4m09.5s 


+54d31m37s 


100.7 


3100 


23200 


25160 


8.8 


(1) 


(2) 


(3) 


(16) 


21 


NGC4670 


12h45ml7.1s 


+27d07m32s 


13.7 


210 


2634 


4470 


8.2 


(1) 


(2) 


(6) 


(12) 


22 


SBS1415+437 


14hl7m01.4s 


+43d30m05s 


4.3 


19.6 






7.6 


(14) 


(5) 


(9) 


23 


Mrk475 


14h39m05.4s 


+36d48m21s 


<2.7 


10.8 






7.9 


(14) 


(5) 




(9) 


24 


CG0563 


14h52m05.7s 


+38dl0m59s 


6.2 


97.0 


870 


1900 


8.7 


(1) 




(6) 




25 


CG0752 


15h31m21.3s 


+47d01m24s 


4.8 


138.9 


837 


1059 




(1) 


(5) 


(6) 




26 


SBS1533+574 


15h34ml3.8s 


+57dl7m06s 


4.2 


53.4 


257 


405 


8.1 


(14) 


(5) 


(6) 


(9) 


27 


Mrkl499 


16h35m21.1s 


+52dl2m53s 


1.5 


33.8 


256 


617 


8.1 


(20) 


(5) 


(6) 


(21) 


28 


Mrk930 


23h31m58.3s 


+28d56m50s 


12.2 


170 


1245 


<2154 


8.1 


(1) 


(2) 


(6) 


(7) 



Ref erenc es. — (1) Condon ct al. (1998), (2) Engclbracht ot al. (2005'). f3l [ Sanders et al.l <2003D. ( 41 [Bergvall fc OstlinI 1I200I1. fSl IWu et al.l 
((20060, (6) 'Moshir & ct al. (1990), (7) Gusova ct aL (2000), (8) Hunt ct al (!2o6j), (9) |Izotov fc ThuanP (1199^. (10) 'Kobulnickv &: Skillmanl (119971)7 
(11) Kobulnickv & .lohnson (1999), (12) Hcckman ct al. (1998), (13) Hunt ct al. (2005a), (14) Hopkins ct al. (2002), (15) Thuan & Izotov (2005), (16) 
[Storchi-Bergmann ct al.i (11994) , ('17) ,Kunth &: Jou bcrt (1985) , (18) .Lerov ct al.i (i2005i) , (19) ilzotov k. Tlmaa (,1998,') , (20) FIRST catalog, (21) .Shi ct al.l 
( 120051) 

the photometry of our sources in Table [T] The uncer- 
tainty in the 22 /im or 24 /im photometry is less than 
5%. The IRAS 60 /im and 100 /im fluxes typically have 
less than a few percent of error, but could go up to ^^15% 
for some of the fainter sources. The rms noise level for 
NVSS is ~0.5mJybeam-i. 

3. RESULTS 

3.1. Mid-IR and FIR to Radio Correlation in BCDs 

The sensitivity and efficiency of the Spitzer Space Tele- 
scope has allowed us to probe the correlation between 
the mid-IR and radio luminosities for a large number of 
galaxies. Using the 24 /<m a nd 70 um MIPS i magin g data 
of the First Look Survey, lAppleton et al.l (|2004[) have 
demonstrated the first direct evidence for the universality 
of the m id-IR/radio and FIR/radio correlation to z^l. 
IWu et al.l (|2009) have also studied the mid-IR/radio cor- 
relation in a sample of star-forming galaxies and found 
that both the 8 /im and 24 /xm luminosity are clearly cor- 
related with the 1.4GIIz radio luminosity. Their sample 
included only a few (3) dwarf galaxies and suggested that 
there may be a slope change for dwarf galaxies, which 
could be due to the lower dust-to-gas ratios and lower 
metallicities of the dwarfs. A detailed analysis of the 
spatial distribution of the infrared to radio correlation 
using Spitzer data on a sample of nearby late type spi ral 
galaxies has been performed bv lMurphv et al.l ()2006ai rbl) 
in which the authors found that the ratio of the mid-IR 
and FIR to radio emission does vary from one region to 




9 10 
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Fig. 1. — The FIR to 1.4 GHz radio luminosity correlation for 
this sample. Our data are presented with filled circles. The num- 
bers next to each symbol correspond to their IDs in Table [l] The 
solid line is the best fit to this sample excluding the upper lim- 
its. For comparison , we have also include the starburst galaxies of 
IBrandl et al.l (|2006l ). indicated with triangles. 

the other. In general though, the dispersion is small and 
the ratio decreases with surface brightness and galacto- 
centric radius. 

Using the data from this sample of low metallicity 
dwarf galaxies we plot in Fig. 1 the radio luminos- 
ity of the sample as a function of the FIR luminos- 
ity. The luminosities of the sources we study span 
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TABLE 2 
Derived Quantities of the Sample 



ID 


Object Name 


Distance^ 


Ll.4GHz 


L24u,m 


Lfir 


q24 


qpiR 






(Mpc) 


(xlO^*^ WHz"^) 


(xIO^Lq) 


(xlO^Lo) 








Uornl 1 

ndi o j_ 1 


oo 


ZuU.U 




I uu.o 


1 


2 46 


2 


IN v_T V_y _L l^U 


21 2 


12 8 




24 4 


11*^ 
1 . j_o 


2 29 


Q 

o 


oJjoUooO-UOZri; 


Oo 


18 7 


C5.0 




9 1 
Z. ID 








9 n 


1 (\ 




1 


u. yo 


2 29 


c; 
O 


llZj w^u 


19 4 


D.U 


Q 9 

y.z 


1 7 


1 «f{ 

i .DD 


Z.oD 


p. 




fi Q 

u. y 


u.u 


U.U 


9 S 
z.o 


1 "^fi 
1 .OU 


z.uo 


7 


He2-10 


12 




9Q Q 
zy .y 




1 77 
J- . M 




Q 

O 


In v^t V_,' Z ( OZ 




91 '^ 


uu.o 


9fin 4 

ZUU.I 


u.oy 


9 nn 

z.uu 


Q 




o. y 


42 5 


u.y 


on Q 
yu. y 


u.uy 


2 34 


10 


IZwl8 


18.2 


0.7 


0.07 




0.48 


1.31*^ 


11 


bi3bUy4U+044 


23 


< 1.5 


0.05 




>0.00 




12 


1NCjOoU77 


3.8 


0.5 


0.9 


3.9 


1.71 


2.89 


13 


Mr K loo 


40.5 


7.9 


1.9 


<7.9 


0.86 


<2.00 


14 


VliZ/w4Uo 


4.3 


0.03 


0.02 


" " " 


1.37 




15 


MrKi4oU 


20.0 


<1.0 


0.8 


<2.1 


>1.38 




16 


T T A T /I /I O 

U M44o 


87.4 


300.6 


188.9 


458.6 


1.23 


2.19 


17 


T T A T /I ^; 1 


13.4 


<0.6 


0.2 




>1.06 




18 


UM462 


13.4 


1.3 


0.8 


2.4 


1.28 


2.29 


19 


SBS1159+545 


50 


<7.0 


0.6 




>0.44 




20 


NGC4194 


41.5 


209.3 


212.0 


586.0 


1.49 


2.45 


21 


NGC4670 


23.2 


8.9 


4.5 


24.3 


1.19 


2.44 


22 


SBS1415+437 


8.7 


0.4 


0.06 




0.66 




23 


Mrk475 


11.2 


<0.4 


0.05 




>0.60 




24 


CG0563 


139 


144.6 


74.4 


320.5 


1.19 


2.35 


25 


CG0752 


90 


47.3 


45.1 


105.3 


1.46 


2.35 


26 


SBS1533+574 


47 


11.3 


4.7 


9.6 


1.10 


1.93 


27 


Mrkl499 


39 


2.7 


2.0 


7.6 


1.35 


2.46 


28 


Mrk930 


77.5 


88.4 


40.5 


<129.0 


1.14 


<2.17 



^ The distances to the galaxies of the sample are adopted from iMoustakas & KennicuttI j2006D when 
available, while the rest of the sources arc calculated from the redshifts taken from NED, assuming a 
ACDM cosmology with Hp ^ 70 km s~^ Mp c~\ H^^O.S and Oa ^ 0.7. For IZwlS, we adopt the 
newly derived distance bv lAloisi et al.l II2007I) . The average uncertainty in the distances is ~5%, mainly 



due to the value of Hq. 
Section 3.4). 



This 



; calculated based on the "equivalent" IRAS 60 and 100 fluxes (see 



nearly 4 orders of magnitudes, but the correlation be- 
tween the FIR and the radio is remarkably tight. The 
scatter in the ratio of the FIR and radio luminosities 
is 0.23 dex, i.e. less than a factor of 2. We per- 
formed a least-squares bisector fit to the data and found: 
log[Li.4GHz(WHz-i)] = l.O9xlog(LFiR/L0)-KlO.75 . For 
comparison, we have inclu ded in the plot the starburst 
galaxies from lBrandl et aLl ([2006) marked with triangles, 
which has an identical slope (within la). This slope of 
1.09±0.07 for the dwarf gala xy data agr ees well with the 
slope of 1.10±0.04 found by[Be3 ([20mh for a sample of 
162 galaxies, as well as the slope of 1.11±0.02 for the 
infrared selected sources from the IR AS Bright Galaxy 
Sample (BGS) ("Cond on et aLl flOQll) . This is also in 
agreement with Hopki ns et al.l ()2002[ ) and would suggest 
that globally our BCDs have a very similar FIR/radio 
correlation to normal galaxies. 

Another way to parameterize the IR/radio correlation 
is to calculate th e ratio of FIR to radio luminosity qfir 
according to the iHelou et al. (igSS.) formula as well as 
the 524 following the definition of lAppleton et al.l ()2004D ^. 
We plot the (724 values for our sample as a function of the 
24/im luminosity of the galaxies in Fig. 2. For this sam- 
ple we find that qfir — 2.4±0.2, consistent with the value 



(2.58S60Mm+SlOOMm)/(3.75 

andSioOAim are in Jy and Fi.4ghz 



^ We use qpiR=log[1.26xlO- 
xlO^^ Fi.4GHz)] whereSeoMT, 

is in Wm-^Hz-i l|Helou et al.l [T985l V We also define 924= 
log (S24um/Si 4 GHz)^ where 824^™, a^nd Si.4ghz are in Jy as in 
lAppleton etall pOM ). 



of nor mal galaxies of qpm = 2.3±0.2 found bv ICondonl 
dHH). When using the mid-IR 24/im data, we find that 
924 =1.3±0.4 (see Table [2]). The standard deviation in 
924 is ^twice that of (/fir- This is not unexpected given 
that the spectrum of star forming galaxies shows sub- 
stantial ly larger variations in spectral slope i n the mid- 
IR (s ee lBrandl et al.l [20061 ) than in the FIR (iDale et all 
120061) . A small change in the geometry of the emitting 
regions would affect Fj^ (24/im) and thus q24 ratio much 
more than the FIR emissio n and qpiR,. A simila r re- 
sult has also been not iced bv lAppleton et all ()2004[ ) and 
iMurphv et al.l ()2006af ) who found a larger dispersion in 
q24 as compared to qvo and suggested that this is prob- 
ably due to a larger intrinsic dispersion in the IR/radio 
correlation at s horte r wavelengths. Interestingly, the 
lAppleton et"all (|2004f ) values for q24 = 0.84±0.28 or the 
k-corrected q24 of 0.94±0.23 are somewhat smaller than 
our results, though consistent within 2a. One possible 
explanation is that the dust temperature of low luminos- 
ity dwarf galaxies tends to peak at shorter wavelength 
than normal star forming galaxies, which would result in 
an elevated 24 fim based luminosity. An alternative ex- 
planation is that due to the peculiar morphology of dwarf 
galaxies, more electrons esc ape from th e galaxy because 
of the cosmic ray diffusion. iBovle et al.l ((2007. ) have also 
found a high q24=1.39±0.02 for the sources they stud- 



We define q70= log (S70mn/Si . 4GHz) where S70/J 



Si.4GHz are in Jy as in lAppleton et ahl 12004 



and 
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Fig. 2. — The 524 as a function of the the 24 /xm monochromatic 
luminosity L24. The symbols are the same as in Fig. 1. 

ied, and postulated that this may be due to a change in 
the mean q24 ratio for objects with Fy(24^m) <lmJy, 
however, their sources are much further away and may 
not be comparable to the low metallicity star forming 
galaxies we study in the local universe. 

3.2. Metallicity and Dust Temperature Effects on qm 

A number of physical parameters, such as the metal- 
licity, dust grain size distribution, and temperature may 
affect the shape of the infrared SED. Consequently, we 
search for correlations between the q ratios of the galaxies 
in our sample with those parameters. Our sample cov- 
ers a metallicity range of 7.2<12-|-log(0/H)<8.9. In low 
metallicity galaxies, the dust content is usually low, even 
thoug h there are notable exceptions such as SBS0335- 
052E (jHouck et al.ll2004bD . thus the emission of UV hght 
might not be fully reprocessed by the dust and re-emitted 
at the infrared wavelengths. However, these galaxies 
might have a quenched synchrotron radiation. This can 
be attributed to various reasons including a lack of super- 
novae remnants which accelate particles producing radio 
emission, the escape of fast cosmic rays from the galaxy, 
etc. These two competing factors, d ust and rad io emis- 
sion, counter balance each other (see lBellll2003f ). BCDs 
though, are typically small, less than ~1 kpc in size, and 
have fewer HII regions with massive stars (w.r.t. nor- 
mal galaxies) which heat the dust and go supernovae 
to produce cosmic rays. As a result the averaging ef- 
fects in the sampling of the properties of the interstellar 
medium which result from the limited spatial resolution 
in dense galactic disks is not so prominent in BCDs. Fur- 
thermore, the small number of statistics may also con- 
tribute to the observed higher dispersion in q24 ratios in 
metal-poor dwarf galaxies. We investigate how metallic- 
ity affects the q24 ratios by dividing the dwarf galaxies 
in our sample into two groups: a lower metallicity group 
with 12-|-log(O/H)<8.0 and a higher metallicity group 
with 12-Hog(O/ H)>8.0. This metallicity threshold was 
selected following iRosenberg et all ()2007[ ) who noticed a 
change in the properties of the star-forming dwarf galax- 
ies they studied around this metallicity value. We find 
a mean q24 =1.1±0.1^^ for the first group and q24 = 

Here we use the standard error of the mean (SEM) to quantify 
the dispersion in the mean value. 



Fig. 3. — The q24 ratio plotted as a function of the oxygen 
abundance of the BCDs. The mean q24 value for all the sources 
is indicated by the solid line. A fit to the low metallicity sources 
(12+log(O/H)<8.0) is indicated by the dashed line; while the dot- 
ted line is the same fit excluding SBS0335-052E. 

1.3±0.1 for the second one. If we calculate the q24 for the 
low metallicity group without including SBS0335-052E, 
we find a much lower q24 = 0.9±0.1. We also notice that 
it appears q24 generally decreases with reduced metal- 
licity for sources with 12-l-log(O/H)<8.0, with SBS035- 
052E as a clear outlier, and the slope flattens out at 
12-Mog(O/H)>8.0 (see Fig. 3). This scatter in q24 at 
higher metallicities could be attributed to the more dis- 
persion in the mid-IR emission in these sources and is 
consistent with a number of other studies that have found 
increased dispersion in the q ratios for very high lumi- 
nosity galaxies (whi c h usually have higher meta llicities) 
(iCondon e t Zl[l99ll: 1^ et al.ll2001t iBelillM)! . 

Hummel et al.l (|1988| ) have found evidence that, 
for a given galaxy, there is a small decrease in 
F(100/im)/S(20cm) when the dust temperature in- 
creases. The latter can be traced by t he ratio of 
F(60/im)/F(100/xm). iRoussel et all (|2003D reached a 
similar result for the starburst galaxies they stud- 
ied, in which they detected an anti-correlation of 
qpiR with F(60/im)/F(100/im). However, perform- 
ing the same test for the q24 ratios on our dwarf 
galaxy sample which spans a range of flux ratios -0.4 
<log(F(60^m)/F(100/;im))<0.2 we find no such clear 
trend. It is conceivable that this is partly because our 
sample is too small to identify a trend due to the intrinsic 
scatter in q24. 

3.3. Star Formation Rate Estimates 

Both the radio and infrared emission can be used to es- 
timate the star formation rates (SFRs) in galaxies. How- 
ever, these correlations depend on a number of parame- 
ters including dust content, optical depth and metallicity. 
This topic and potential caveat s have been di scussed ex- 
tensi vely in the literature fsee ICondon|[l992t iKennicuttj 
Il998l and references therein). For dwa rf galaxies, the 
topic has been addressed in detail by I Hopkins et al.l 
(|2002| ). The recent wealth of data from Spitzer has also 
provided sufficient motivation to establish a calibr ation 
for the SFR using the infra r ed. IWu et al.l ^005) and 
more recently Calzetti et al.l (|2007t ) have explored this 
topic using a large sample of star forming regions and 
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nearby galaxies. 

In Fig. 4 we plot the IR estimate d SFR for our sa mple 
using the calibration s proposed by IWu et all (|2005( ) and 
ICalzetti et al.l ()2007D respectively: 



SFR24(A/oyr- 



6.66 X lOSLfTi 



10.885 



(1) 



(2) 



SFR24(M02/r-i) = 1.27 x IQ-^'^lL^Aiergs-^)] 

as a function o f the well-known radio to SFR formula 
of lCondonI (Il993 ): 



SFRi.4GHz = 5.5 X .'aimm n 

4.6 X 10^^(WHz^^) 

The IWu et all (|2005l ) work is based on a global cor- 
relation without separating the low metallicity sources 
from metal rich galaxies. As can be seen in Fig. 4, 
a good agreement exists between the radio and IR es- 
timated SFRs (indicated by diamonds) . If w e use the 
more recent calibration bv lCalzetti et al.l ()2007 ^ on SFRs 
from 24 /im luminosities, we find that most of the mid- 
IR estimated SFRs (marked as filled circles) are located 
below the 1:1 proportionality line, and thus they are con- 
sistently lowe r than both the SFRs estimated from the 
radio or from IWu et"aI1 (|2005f ). This is not unexpected 
since eq. 2 is calibrated based on the high metallicity 
sources that have significant 24 /im emission, while most 
of our so urces are meta l -poor galaxies. We should also 
note that ICalzetti et all (|2007f ) measure the SFR in indi- 
vidual Hii regions in apertures within disks and subtract 
a "disk background" , which could partly explain the de- 
viations we see from our analysis. These authors have 
also mentioned that the SFRs of low metallicity galax- 
ies would be underestimated by a factor of 2-4 depend- 
ing on how metal-poor the galaxies are. This deviation 
from a linear correlation is likely d ue to the lower opac - 
ities for decreasing metal content () Walter et al.l l2007f ) . 
When we fit the filled circles on Fig. 4, we find 
that log(SFR[IR])=0.94xlog(SFR[1.4GHz]/3.98), which 
is consist ent with the m e tallic ity correction factor sug- 
gested by ICalzetti et al.l (I2007D. We should stress once 
more though that as "'Bell' (2003) has noted, the good 
agreement between the IR and radio estimates of SFRs 
does not necessarily mean that these are the "true" SFRs 
for the dwarf galaxies we study, but rather the com- 
petition of the lower dust content and suppressed syn- 
chrotron emission balances each other. 

3.4. The Two Extremes: IZwl8 and SBS0335-052E 

Despite this overall agreement between the IR/radio 
correlation in our dwarf sample and the correspond- 
ing values of normal galaxies, the two lowest metallic- 
ity galaxies in our sample, IZwl8 and SBS0335-052E, 
deviate markedly from the correlations and in opposite 
directions (see Fig. 3"'^^ ). As the two most well-studied 
BCDs, it is interesting to inspect the IR/radio properties 
of these two galaxies and a comparison of their physical 
parameters can be found in Table [3l 

12 In Fig. 2 the "outlier" status of SBS0335-052E does not ap- 
pear as prominent as in Fig. 3, because its high 24 /xm luminosity 
places it in a region of the plot where other sources with higher 
metallicity and higher luminosity are located. 





a. 

^ -1 



'I I I I ' 



'I I ' 



Wu et ol. (2005)0 
Colzett: et ol. (2007) • 



-4-3-2-10 1 2 3 
log (SFR [1.4GHz])(Mo/yr) 

Fig. 4. — A plot of the SFR based on the 24 fim luminosity as 
a function of the corresponding SFR estimated from the 1.4 GHz 
radio continuum emission for each galaxy in our sample. The open 
diamo nds indicat e the SFRs calculated using the 24 fj,m fluxes fol- 
lowi ngjWiret al.l EOOS.') (e a.l). the filled circles are those using eq.2 
fromCalzct ti et al.l II2007I) . The solid line is the 1:1 proportionality 
line. The dotted line is a fit to the filled circles. 



TABLE 3 

Comparison of SBS 0335-052E and I Zw 18 





SBS 0335-052E 


IZwl8 


Ratio^ 


124-log(0/H) 


7.32 


7.17 


1.4 


D (Mpc) 


58 


18 


3 


R (pc) 


560*= 


870^= 


~0.6 


SFR(M0yr-i) 


7x10^1 


5x10-2 


14 


SNR(yr-i) 


6x10-3 


2x10-4 


30 




5.6x10'''* 


1.6x10^ 


35 


L24fim(LQ) 


9x10^ 


7x10*^ 


120 


Lir(Lq) 


^4x10^^ 


~4xl0'' 


100 


Ll.4GHz(WHz-l) 


2x10^0 


8x10^'' 


2.5 


M(Ks)(mag) 


-18.3 


-16.1 


~80<= 



^ Ratio of the relative quantities in SBS0335-052E to 
IZwlS.^ This is the size of the six super star clusters in 
SBS0335-052E, which has a size of 2" jThuan et al.||199iD . 
At a distance of 58 Mpc, 1 ~280 pc.^ This indicate the 
largest projected linear extend of the total star- forming re- 
gions m this galaxy (•'^10". lHunt et al.|[2005bll . At a distance 
of 18 Mpc, l" ~ 87 pc.** Flux in l" slit multiplied by a factor 
of 2 to correct for aperture effects.^ This is the approximate 
mass ratio of the two sources assuming that their mass scales 
with K-band luminosity and that both sources have the same 
M/Lkb fsee lBell fc de JondHoOlf) . 



As indic ated by the Spit zer IRS spectrum of SBS0335- 
052E ( Houck et al.ll2004bf ). the SED of the galaxy has 
an unusual l y flat m id-IR slope that peaks at ~28Mm. 
Hunt et all (|2005af) used DUSTY models (jlvezic et al.l 



Tool " to fit its SED and found a qFiR value of 2.0 for 
this galaxy, lower than the typical qFiR of '~2.3 expected 
for late type systems. These calculations are based on a 
1.4GIIz continuum flux of '-~^0.46mJy, which c orresponds 
to th e compact part (6 ") of SBS0335-052E (jPale et all 
1200 ll : iHunt et al.ll2004[ ). Similarly, if we simply calcu- 
late the q ratio independent of any modelling work and 
only use the observable parameters, we find a qvo = 
2.1±0.1 (assui ning that the MIPS 70 u m flux density is 
51.1±4.8mJv. lEngelbracht et al.|[2007n . This is in good 
agreement with the qyo = 2.16±0.17 for the sources in 
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the First Look Survey studied bv lAppleton et al.l (|2004D . 

However, when we examine the q24 for SBS0335-052E, 
adopting a MIPS 24 /um flux of 66 mJy, we find a q24 ra- 
tio of 2.2±0.1. We see a 2(T excess in the q parameter 
compared with the mean q24 of 1.3±0.4 in the dwarf sam- 
ple. If this excess is real, it would be consistent with the 
unique shape of the SED of the galaxy, which suggests 
that the dust grain distribution is dominated by small 
grains with temper atures of ~150K an d a total dust 
mass of - 10^ MrT. (|Houck et al.ll2004bD . Whether the 
large grains were never created or they were destroyed 
by shocks is unknown. The elevated q24 r atio could also 
be due to th e lack of synchrotron emission. iRoussel et al.l 
(|2003l |2006[ ) have shown their study of three starburst 
galaxies, all of which have qpiR more than 3cr higher 
the mean qpiR and categorized them to be nascent star- 
bursts. SBS0335-052E though, could also be a candidate 
of this group. It is based partly on the high q24 value, as 
well as its low Ha luminosity as compared to the infrared 
or radio luminosities (see Table [3]) . Finally, 'Hu nt etldl 
(|2b04 ) have shown the existence of free-free absorption 
in the radio spectrum of SBS0335-052E, which is usually 
caused by young, dense and heavily embedded clusters^^. 

IZwl8, despite its similar metallicity to SBS0335- 
052E (12H-log(0/H)=7.17 and 7.32 respectively), shows 
a rather different q ratio. It also has an IR luminosity of 
LiR~lO^L0, just 1% of the Ljr in SBS0335-052E. In low 
luminosity galaxies, radio emission is known to dec rease 
faster than dust emission (jPevereux fc EalesllT989h . and 
high IR/radio ratios have been observed in low luminos- 
ity dwarf galaxies (,Klein et a l. 1984). In IZwl8 though, 
we find a rather low infrared to radio ratio. The q24 
is found to be 0.5±0.1, nearly 2(t lower than the aver- 
age value of our sample. At longer wavelengths, IZwl8 
is faint and no IRAS FIR dat a are available. However, 
as noted bv IWu et all (|2007af ) (see their Fig. 5), it has 
a very similar 5-38/im cont inuum slope to the typical 
starburst galaxy NGC7714 (|Brandl et al.l l2004f) . If we 
were to assume that this similarity extends to the FIR 
and scale down the IRAS 60 and 100 /im flux densities 
of NGC7714 by a factor of 375 so that its correspond- 
ing 22 /im flux density matches the one of IZwl8, we 
flnd its "equivalent" /i?^ 5 60 and 100 /im flux density to 
be 29.8 and 32.8 mJy respectively. This would result in 
qpiR, of 1.3 for the galaxy, which deviate s from the av er- 
age qpiR ratio of 2.3±0.2 calibrated by Co ndon! ()1992( ) by 
nearly 5(t. If we w ere to use the MIPS 70 /im detection of 
34mJy for IZwl8 (jEngelbracht et al.ll2007,) and calculate 
its q ratio, we would flnd q7 n=1.4, again 5a away from 
the average qvo sugg ested bv lAppleton et al.l (|2004D . 

A number of plausible scenarios were considered to ex- 
plain this r esult, though non e appears convincing. As 
discussed in lWu et"al] ()2007aD . because the optical depth 
of IZwl8 is small, it could be that a significant fraction 
of the UV light has leaked out without being absorbed 
by the dust, thus resulting in a damped 24 /im emission. 
Alternatively, IZwl8 may simply have an unusually high 
radio luminosity at 1.4 GHz. Could it be that IZwl8 is 
observed at a special moment right after the explosion 

Note, that if we could properly account for the self-absorption 
of the radio continuum, the intrinsic q ratio would decrease. How- 
ever, the exact fraction for absorption is not known and more data 
at u <1.5GHz would be needed to determine that number. 



of a supernovae event? Radio su pernovae fade by more 
than a factor of 10 within -3yr (|Chevalieilll982l) . The 
radio continuum obser vations of IZwl8 s pan over a pe- 
riod o f more than 5yrs ()Hunt et a l. 2005b: ICannon et al.l 
|2005| ) but show no varia tion, thus this i s not likely. 

The recent results by iMurphv et al.l p006b ) provide 
another possible scenario. These authors analyzed a sam- 
ple of nearby spiral galaxies and found that systems with 
higher disk-averaged SFRs (S sfr) have usually experi- 
enced a recent episode of enhanced star formation. As a 
result they contain a higher fraction of young cosmic ray 
electrons that have traveled only a few hundred parsecs 
from their acceleration sites in supernova remnants. This 
is perhaps the case for IZwl8. However, other sources 
in our sample also have elevated E sfr and show no 
radio excess. It could also be that the extremely low 
metal abundance of IZwl8, 0.15 dex lower than that of 
SBS0335-052E, is below a critical threshold. The newly 
formed stars in the unpolluted interstellar medium may 
produce, and subsequently heat, less dust than electrons 
which are accelerated and contribute in the radio emis- 
sion. We also note that the rate of SNR in IZwl8 is a 
factor of ~30 lower than that in SBS0335-052E (see Ta- 
ble O) but its morphology is much mor e disturbed with 
filamentary structure and outflows (see llzotov fc Thuaiil 
2004). If we were to assume that both sources have sim- 
ilar mass-to-light ratios and estimate their mass from 
their K-band luminosities, we find that the mass of IZwl8 
is almost a factor of 80 less than that of SBS0335-052E. 
Thus the rate of SNR normalized with mass in IZwl8 
would be ~2.5 times that of SBS0335-052E. Could it be 
that the disturbed morphology of IZwl8 in addition to 
its low metal content that it created the conditions in the 
interstellar medium for this abnormally high radio flux? 
The question remains open. 

4. CONCLUSIONS 

We have studied the mid-IR and FIR to radio correla- 
tion in a sample of dwarf star-forming galaxies spanning 
a metallicity range from 7.2<12-f log(0/H)<8.9, using 
Spitzer IRS/MIPS, as wefl as radio 1.4 GHz data ob- 
tained from the literature. The BCD sample appears 
to follow the same FIR/radio correlation as normal star 
forming galaxies. The analysis based on mid-IR and ra- 
dio data reveals a similar correlation, though the scatter 
is larger, probably due to an intrinsically higher varia- 
tion in the 15-30/im SEDs of dwarf galaxies. When com- 
paring the q24 ratios with metallicity or effective dust 
temperature, we flnd no strong correlation, though there 
is a general trend of lower q ratios at lower metallicity 
for galaxies with 12-|-log(O/H)<8.0 and the correlation 
flattens out toward higher metallicity. In general the 
SFRs estimated from the radio 1.4 GHz continuum and 
the mid-IR data are in good agreement. Two extremely 
metal poor BCDs, IZwl8 and SBS0335-052E appear to 
deviate from the average q24 by ~2cr, with one displaying 
a radio excess and the other an infrared excess. 
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